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ROOT'S BARRIER: CONSTRUCTION, OPTIMALITY AND 
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Recent work of Dupire and Carr and Lee has highlighted the im- 
portance of understanding the Skorokhod embedding originally pro- 
posed by Root for the model-independent hedging of variance options. 
Root's work shows that there exists a barrier from which one may de- 
fine a stopping time which solves the Skorokhod embedding problem. 
This construction has the remarkable property, proved by Rost, that 
it minimizes the variance of the stopping time among all solutions. 

In this work, we prove a characterization of Root's barrier in terms 
of the solution to a variational inequality, and we give an alternative 
proof of the optimality property which has an important consequence 
for the construction of subhedging strategies in the financial context. 

1. Introduction. In this paper, we analyze the solution to the Skorokhod 
embedding problem originally given by Root [33] , and generalized by Rost [35] . 
Our motivation for this is recent work connecting the solution to this prob- 
lem to questions arising in mathematical finance — specifically model-indepen- 
dent bounds for variance options — which has been observed by Dupire [16], 
Carr and Lee [5] and Hobson [19]. The financial motivation can be described 
as follows: consider a (discounted) asset which has dynamics under the risk- 
neutral measure 

where the process at is not necessarily known. We are interested in vari- 
ance options, which are contracts where the payoff depends on the realized 
quadratic variation of the log-price process: specifically, we have 

d{ln St) = at dWt-la^t 
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and therefore 



■T 



{\nS)T 



a^dt. 



An option on variance is then an option with payoff F{{lnS)T)- Important 
examples include variance swaps, which pay the holder {In S)t — K, and 
variance calls which pay the holder ((IniS')^ — We shall be particularly 
interested in the case of a variance call, but our results will extend to a wider 
class of payoffs. Let dXt = Xt dWt for a suitable Brownian motion Wt and 
we can find a (continuous) time change Tf such that St = X-j-^ , and so 



Now suppose that we know the prices of call options on St with maturity T, 
and at all strikes (recall that at is not assumed known). Then we can derive 
the law of St under the risk-neutral measure from the Breeden-Litzenberger 
formula. Call this law This suggests that the problem of finding a lower 
bound on the price of a variance call (for an unknown at) is equivalent to 

(1.1) find a stopping time r to minimize E(r — -fC)+, subject to C{Xt-) = fJ-. 

This is essentially the problem for which Rost has shown that the solution 
is given by Root's barrier. [In fact, the result trivially extends to payoffs of 
the form F{{lnS)T) where F(-) is a convex, increasing function.] 

In this work, our aim is twofold: first, to provide a proof that Root's 
barrier can be found as the solution to a particular variational inequality, 
which can be thought of as the generalization of an obstacle problem; second, 
we show that the lower bound which is implied by Rost's result can be 
enforced through a suitable hedging strategy, which will give an arbitrage 
whenever the price of a variance call trades below the given lower bound. 
To accomplish this second part of the paper, we will give a novel proof of 
the optimality of Root's construction, and from this construction we will be 
able to derive a suitable hedging strategy. 

The use of Skorokhod embedding techniques to solve model-independent 
(or robust) hedging problems in finance can be traced back to Hobson [18]. 
More recent results in this direction include Cox, Hobson and Obloj [10], 
Cox and Obloj [11] and Cox and Obloj [12]. For a comprehensive survey 
of the literature on the Skorokhod embedding problem, we refer the reader 
to Obloj [27]. In addition, Hobson [19] surveys the literature on the Sko- 
rokhod embedding problem with a specific emphasis on the applications in 
mathematical finance. 



drt 




dt. 
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Variance options have been a topic of much interest in recent years, both 
from the industrial point of view, where innovations such as the VIX index 
have contributed to a large growth in products which are directly depen- 
dent on quantities derived from the quadratic variation, and also on the 
academic side, with a number of interesting contributions in the literature. 
The academic results go back to work of Dupire [15] and Neuberger [26], 
who noted that a variance swap — that is, a contract which pays (ln5')T, can 
be replicated model-independently using a contract paying the logarithm of 
the asset at maturity through the identity (from Ito's lemma) 



More recently, work on options and swaps on volatility and variance, (in 
a model-based setting) includes Howison, Rafailidis and Rasmussen [20], 
Broadie and Jain [3] and Kallsen, Muhle-Karbe and Voss [21]. Other work 
[22, 23] has considered the differences between the theoretical payoff ((ln5)r) 
and the discrete approximation which is usually specified in the contract 
[}2k^^i'^{k+i)s/Sk5)'^]- Finally, several papers have considered variants on 
the mo del- independent problems [5, 6, 13] or problems where the modeling 
assumptions are fairly weak. This latter framework is of particular inter- 
est for options on variance, since the markets for such products are still 
fairly young, and so making strong modelling assumptions might not be as 
strongly justified as it could be in a well-established market. 

The rest of this paper is structured as follows: in Section 2 we review 
some known results and properties concerning Root's barrier. In Section 3, 
we establish a connection between Root's solution and an obstacle problem, 
and then in Section 4 we show that by considering an obstacle problem in 
a more general analytic sense (as a variational inequality), we are able to 
prove the equivalence between Root's problem and the solution to a vari- 
ational inequality. In Section 5, we give a new proof of the optimality of 
Root's solution and in Section 6 we show how this proof allows us to con- 
struct model-independent subhedges to give bounds on the price of variance 
options. 

2. Features of Root's solution. Our interest is in Root's solution to the 
Skorokhod embedding problem. Simply stated, for a process {Xt)t>o, the 
Skorokhod embedding problem is to find a stopping time r such that Xr ~ fi. 
In this paper, we will consider first the case where Xq = 0, and Xt is a contin- 
uous martingale and a time-homogeneous diffusion, and later the case where 
Xq ~ z/, is a centred, square integrable measure. In such circumstances, it 
is natural to restrict to the set of stopping times for which {XtAT)t>o is a 
uniformly integrable (UI) process. We will occasionally call stopping times 
for which this is true UI stopping times. In the case where /i is centered and 



(1.2) 
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has a second moment and the underlying process X is a Brownian motion 
(or more generally, a diffusion and martingale with diffusion coefficient a 
such that fJ^ > e for some strictly positive constant e), this is equivalent to 
the fact that Er < oo. For the case of a general starting measure, there is a 
natural restriction on the measures involved, which is that we require 

(2.1) oo > U,,(2;) := - / \y - x\v{dy)>- / |y - x|/i(dy) =: U^(x), 

for all X S M. This assumption implies that m := J xi'{dx) = J xfi{dx); see 
Chacon [7]. By Jensen's inequality, such a constraint is clearly necessary for 
the existence of a suitable pair and /i; further, by Rost [34], it is the only 
additional constraint on the measures we will need to impose. We shall write 

(2.2) S{f-i) = {r:r is a stopping time,^,- ~ fi, {Xt/\T)t>o is UI}. 

There are a number of important papers concerning the construction of 
Root's barrier. The first work to consider the problem is Root [33], and 
this paper proved the existence of a certain Skorokhod embedding when Xt 
is a Brownian motion. Specifically, Root showed that if Xt is a Brownian 
motion with Xq = 0, and /i is the law of a centered random variable with 
finite variance, then there exists a stopping time r, which is the first hitting 
time of a barrier, which is defined as follows: 

Definition 2.1 (Root's barrier). A closed subset B of [— oo,+oo] x 
[0, +oo] is a barrier if: 

(1) (x, +oo) e i? for all X G [— oo, +oo]; 

(2) {±oo,t) e B for ah t G [0,oo]; 

(3) if {x,t) G B, then (x, s) € B whenever s> t. 

We provide representative examples of barrier functions in Figure 1. 

In a subsequent paper Loynes [24] proved a number of results relating 
to barriers. From our perspective, the most important are, first, that the 
barrier B can be written as B = {{x,t) :t > R{x)}, where ii:]R — )■ [0,oo] is a 
lower semi-continuous function (with the obvious extensions to the definition 
to cover R{x) = oo); we will make frequent use of this representation. In 
addition, Loynes [24], Theorem 1, says that Root's solution is essentially 
unique: if there are two barriers which embed the same distribution with a 
UI stopping time, then their corresponding stopping times are equal with 
probability one. The case where two different barriers can occur are then only 
the cases where, say R{xo) = for xq > 0, and then R{x) is undetermined 
for all x> xq. 

The other important reference for our purposes is Rost [35]. This work 
vastly extends the generality of the results of Root and Loynes, and uses 
mostly potential-theoretic techniques. Rost works in the generality of a 
Markov process Xt on a compact metric space E, which satisfies the strong 
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Fig. 1. Examples of Root's barriers: the representation (a) is an example of a "nice" 
barrier, where some explicit calculations can be made (see Example 5.6); in (b) we observe 
some of the nastier features which a harrier may possess, including spikes, corresponding 
to atoms of the distribution and regions in which the harrier can be unbounded. 

Markov property and is right-continuous. Then Rost recalls (from an orig- 
inal definition of Dinges [14] in the discrete setting) the notion of minimal 
residual expectation: 

Definition 2.2. A stopping time r* G S{ii) is of minimal residual ex- 
pectation if, for each t S , it minimizes the quantity 



over all r G S{^). 

Then Rost proves that [under (2.1)] there exists a stopping time of mini- 
mal residual expectation [35], Theorem 1, and that the hitting time of any 
barrier is of minimal residual expectation [35], Theorem 2. Finally, Rost also 
shows that the barrier stopping times are, to a degree, unique [35], Corol- 
lary to Theorem 2. The relevant result for our purposes (where there is a 
stronger form of uniqueness) is the corollary to Theorem 3 therein, which 
says that if Xt is a process for which the one-point sets are regular, then 
any stopping time of minimal residual expectation is Root's stopping time. 
The class of processes for which the one-point sets are regular include the 
class of time-homogenous diffusions we consider. 

Note that a stopping time is of minimal residual expectation if and only if, 
for every convex, increasing function F[t) (where, without loss of generality, 
we take -F(O) = -P+(0) =0), it minimizes the quantity 




X 



(a) 



(b) 





this fact beini 



■g 



a consequence of the above representation. 
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There are a number of important properties that the Root barrier pos- 
sesses. First, we note that, as a consequence of the fact that B is closed 
and the third property of Definition 2.1, the barrier is regular (i.e., if we 
start at a point in the barrier, we will almost surely return to the barrier 
instantly) for the class of processes we will consider (time-homogeneous dif- 
fusions) this will have important analytical benefits. Second, for a point 
{x,t) ^ B, we know that if the stopped process at time t is at x, then we 
have not yet reached the stopping time for the embedding. This will help in 
our characterization of the law of the stopped process (Lemma 3.2). 

In the rest of this paper, we will then say that a barrier is either a lower 
semi-continuous function ii:M — )• M+, with R{0) ^ 0, or the complement 
of the corresponding connected open set D = {{x,t) :0 < t < R{x)} = R x 
(0,oo) \ B. As noted above, by Loynes [24] this is equivalent to the barrier 
as defined in Definition 2.1. We will define the hitting time of the barrier as: 
tj:) = inf{t > : {Xt,t) ^ D}. Note that the barrier B is closed and regular, so 
that (Xro,TD) G B and P(^''*)(tz) = 0) = 1 whenever {x,t) G B, where P^^''*) 
is the law of our diffusion started at x at time t. 

Finally, we give some examples where the barrier function can be explicitly 
calculated. We note that explicit examples appear to be the exception, and 
in general are hard to compute. First, if ^ is a Normal distribution, we easily 
see that R{x) is a constant. Second, if fi consists of two atoms (weighted 
appropriately) at a < < 6 say, the corresponding barrier is 



In this example, observe that the function R{x) is not unique: we can choose 
any behavior outside [a,b], and achieve the same stopping time. Second, 
we note the that there are even more general solutions to the Skorokhod 
embedding problem (without the uniform integrability condition) since there 
are also barriers of the form 



which will embed the same law (provided ta,tb > are chosen suitably), 
but which do not satisfy the uniform integrability condition. In general, a 
barrier can exhibit some fairly nasty features: consider, for example, the 
canonical measure on a middle third Cantor set C (scaled so that it is on 
[—1,1]). Root's result tells us that there exists a barrier which embeds this 
distribution, and clearly the resulting barrier function must be finite only 
on the Cantor set; however, the target distribution has no atoms, so that 
the "spikes" in the barrier function can not be isolated (i.e., we must have 
liminfy-|-2. R{y) = liminfy^a, R{y) = R{x) for all x E (— 1, 1) n C). 




oo 



X ^ (a, 6), 
X G (a, b). 




X = b, 
X ^ {a, b} 
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3. Connecting Root's problem and an obstacle problem. We now con- 
sider alternative methods for describing Root's barrier. We will, in general, 
be interested in this question when our underlying process Xf is a solution 
to 

(3.1) dXt = a{Xt)dWu Xo^u, 

for a Brownian motion {Wt)t>Q, and we will introduce our concepts in this 
general context. Initially, we assume that cr:M— )-]R satisfies, for some posi- 
tive constant K, 

(3.2) \a{x)-a{y)\<K\x-y\; 

(3.3) 0<(j2(x) < J^(l+x2); 

(3.4) cr is smooth. 

Recall that for the financial application we are interested in, we want the 
specific case a{x) = rc to be included. Clearly, this case is currently excluded; 
however, we will show in Section 4.3 that the results can be extended to 
include this case. 

From standard results on SDEs, (3.2) and (3.3) imply that the unique 
strong solution X"" of (3.1) with v = 5a is a strong Markov process with 
generator ^a^dxx for any initial value a G M. Moreover, (3.4) implies that the 
operator L := ^a^dxx — dt is hypoelliptic; see Stroock [36], Theorem 3.4.1. 

We will write Root's Skorokhod embedding problem as: 

SEP(c7, i^, //): Find a lower-semicontinuous function R{x) such that the 
domain D = {{x,t) :0 < t < R{x)} has X-rjy ~ n, and {Xt/\Tjy)t>o is a UI 
process, where ly is the initial law of Xt, and a the diffusion coefficient. 

Our aim is to show that the problem of finding R is essentially equivalent 
to solving an obstacle problem. Assuming that the relevant derivatives exist, 
we shall show that the problem can be stated in the following way: 

OBS(cr,z^,/i): Find a function u{x,t) £ C^'^(IR x E+) such that 



(3.5a) Uj/(x) = ti(x, 0), 

(3.5b) 0>Vf,{x)-u{x,t), 

(3.5c) o>^(x,i)-^cT(x)2|^(x,i), 

(3.5d) ( ^(^'0 - 2^(^)'a^(^'*) ) (Um(^) - ^(^.0) = 0' 



where (3.5c) is interpreted in a distributional sense — that is, we require 
jjy(l){x)^{x,t) + ia(x)2|^(x,t)0'(x)^ dx < 
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whenever (p € is a nonnegative function. Condition (3.5d) can be inter- 
preted more generally as requiring 

du , . 1 . .od'^u . , 

in a distributional sense whenever {\]^[x) — u{x,t)) ^ 0. However, this is an 
open set, and from the hypoellipticity of the operator L, if this holds in a 
weak sense, it will hold in a strong sense. Hence ^^{x, t) would be continuous 
even if we were only to require (3.5d) to hold in a distributional sense. 

In general, we do not expect u to be sufficiently nice that we can easily 
interpret all these statements, and one of the goals of this paper is to give 
a generalization of 0BS((T, z^, /x) that will make sense more widely. Cases in 
which u may not be expected to be C^'^ include the case where fi contains 
atoms (and therefore is not continuously differentiable) . In addition, we 
specify this problem in C^'^ since, in general, we would certainly not expect 
the second derivative to be continuous on the boundary between the two 
types of behavior in (3.5d). 

Theorem 3.1. Suppose D is a solution to SEP{a,i', fi) and is such that 
= -E|XtAro -x\ G C^'^(]R X M+). 
Then u solves OBS{a, v, /i) . 

This gives an initial connection between 0BS((7, fi) and SEP((7, z^, fi). We 
roughly expect solutions to Root's problem to be the unique solutions to the 
obstacle problem (of course, we do not currently know that such solutions 
exist or, when they do, are unique). This suggests that we can attempt 
to solve the obstacle problem to find the solution D to Root's problem. 
In particular, given a solution to 0BS((7, z^, we can now identify D as 
D = {{x,t):\Jf,{x) < u{x,t),t> 0}. 

Lemma 3.2. For any {x,t) G D, F{XtATD G dx) = F{Xt G dx,t< td). 

Proof. By the lower semi-continuity of R, since (x,t) G D, there exists 
/i > such that 

{x - h,x + h) X [Q^t + h) C. D, 

and hence, for any y £ (x — h,x + h), R{y) > t. On the other hand, if to < t, 
we have 

R{XrJ<TD<t, 

and hence, X^-^ ^ {x — h,x + h). Therefore, 

¥{Xt^rn G dx) = ¥{Xt edx,t< td) + r{Xr^ G dx, t > td) 

= ¥{Xtedx,t<TD). □ 
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Lemma 3.3. The measure corresponding to C{Xt]t < T£)) has density 
p^{x,t) with respect to Lebesgue on D, and the density is smooth and satis- 
fies 

This result appears to be standard, but we are unable to find concise 
references. We give a short proof based on [32], Section V.38.5. 

Proof of Lemma 3.3. First note that, as a measure, C{Xt;t < td) is 
dominated by the usual transition measure, so the density p^{x,t) exists. 

Let {xQjto) be a point in D, and we can therefore find an e > such that 
A = {xq — e,XQ + e) x (io — £,to + e) satisfies A<^ D. Then let / be a smooth 
function, supported on A, and by Ito's lemma, 

fiXtAm , t) = f{Xo, 0) + ^* (XsAm . s) dX^ 

Since / is compactly supported, taking t > to + e, the two terms on the 
left disappear, and the first integral term is a martingale. Hence, taking 
expectations, and interchanging the order of differentiation, we get 

/7p"fe..)(^.W^^ + |)/(y.»)dyd. = o. 

Interpreting p^{y,s) as a distribution, we have 

for {x,t) G A, and since the heat operator is hypoelliptic, we conclude that 
p^{x,t) is smooth in A (e.g., Stroock [36], Theorem 3.4.1). □ 

We are now able to prove that any solution to Root's embedding problem 
is a solution to the obstacle problem. 

Proof of Theorem 3.1. We first observe that u(x,0) = — ]E|Xo — x\, 
and Xq ~ i/, so that u{x,0) = — f \y — x\i>{dy) and (3.5a) holds. Second, since 
{Xt/\r^)t>o is a UI process, by (conditional) Jensen's inequality, 

U{x,t) = -E|X - XtArcl > -nn\x - Xr^WJ^tAro]] = ^ , 

and (3.5b) holds. 
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We now consider (3.5c). Suppose (x,t) € D, and note that 

du 

(3.6) — = l-2F{XtAro<x): 

and therefore (in D) by Lemma 3.3 the function u has a smooth second 
derivative in x. Further, we get 

f)2„, ft 

a{x)'^pi:){x, s) ds 



2], 9^^^'^)^^ 



(3.7) 



HmE 

£4,0 



1 rt/\rD 

2e Jo ^(^)^^[^-'^<^''<^+=l^^ 



-E\x - XtArol + k 



where L'f is the local time of the diffusion at x. It follows that u satisfies 
(3.5c) on D, and in fact attains equality there. On the other hand, if {x,t) ^ 
D, it follows from the definition of the barrier that if r/j > t, the diffusion 
cannot cross the line {{x,s) :s>t} in the time interval [t,TD), and hence 

Therefore, for t > R{x), 

E\x - Xt^m \ = \x\+ ELf^,^ = |x| + EL^^ = E\x - X^J, 

where the last equality holds because td is a UI stopping time. So (3.5b) 
holds with equality when (x,t) ^ D. In particular, we can deduce that either 
(if {x,t) € D) we have equality in (3.5c), or we have equality in (3.5b), in 
which case (3.5d) must hold. It remains to show that (3.5c) holds when 
{x,t) ^ D. However, to see this, consider {x,t) ^ D, and note first that 
u{x, s) = u{x, t) = U^(x) whenever s>t, since (x, s) ^ D. Hence ^{x, t) = 0. 
It is straightforward to check that u{x,t) is concave in x, and therefore that 
0(x,t) < 0, and (3.5c) also holds. □ 

This result connects Root's problem and the obstacle problem under a 
smoothness assumption on the function u. However, ideally we want a one- 
to-one correspondence. We know from the results of Rost [35] that there 
always exists a solution to SEP(cj, z^, ;u), and from Loynes [24] that the solu- 
tion is unique. Our aim is to show that a similar combination of existence 
and uniqueness hold for the corresponding analytic formulation. As already 
noted, we cannot make a strong smoothness assumption on the function 
u{x,t) as required by OBS(cj, z^, ;u), and so we need a weaker formulation of 
this problem. Generalizations of the obstacle problem are well understood, 
and commonly called variational inequalities. In the next section, we will re- 
formulate the obstacle problem as a variational inequality, and we are able 
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to state a problem for which existence and uniqueness are known due to 
existing results. 

4. Root's barrier and variational inequalities. We now study the relation 
between Root's Skorokhod embedding problem and a variational inequality. 
Our notation and definitions, and some of the key results which we will use, 
come from Bensoussan and Lions [1]. 

4.1. Variational inequalities. We begin with some necessary notation 
and results concerning evolutionary variational inequalities. Given a con- 
stant A > and a finite time T > 0, we define the Banach spaces H"^'"^ C 
L2(M) and L'^{0,T; H"^'^) with the norms 



\w\ 



\g\\jjm.>. 



L2(0,T;H'".^) 



E 

k=0 



'2\\x\ 



(x) 



dx; 







where the derivatives Q^i^) are to be interpreted as weak derivatives — that 
is, Q^{x) is defined by the requirement that 



Qk 



dx^ 



for all (j) G C^(M), and is the set of compactly supported, smooth func- 
tions on M. In particular, the spaces H"^''^ and L'^{0,T; H"^'^) are Hilbert 
spaces with respect to the obvious inner products. In addition, elements of 
the set H^'"^ can always be taken to be continuous, and is dense in H^'^] 
see, for example, Friedman [17], Theorem 5.5.20. 

For functions a(x,t),b{x,t) G L°°{M. x (0,T)), we define an operator 



ax{t;v,w) 



-2X\x\ 



dv dw dv 
a{x,t) — — + b{x,t)—w 
ox Ox ox 



dx, 



for v,w ^ L'^{0,T; H^''^). Moreover if da/dx exits, we define, for v ^ H 



2,A 



d 



dv\ 



dv 



A{t)v = -—[a{x,t)— j +(6(x,t) + 2Ao(x,t)sgn(x)) — . 
And finally, for v,w ^ H^''^, 

{v,w)x -- 

so that, for suitably differentiable test functions <^(x) and v G i/^''^, 

{(j),A{t)v)x = ax{t]v,(t)). 



e '^^^^^vwdx, 
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Then we have the foUowing restatement of Bensoussan and Lions [1], The- 
orem 2.2, and Section 2.15, Chapter 3: 

Theorem 4.1. For any given A > and T > 0, suppose: 

(1) a,b,^ are bounded on M x (0,T) with a{x,t) >a a.e. in M. x (0,r) 
for some a > 0; 

(2) V, f G L^O, T; H^''), v G H^'\ v>m ; 

(3) the set 

w{t) > ip{t) a.e. t in [0,T] 

is nonempty, where [H^''^)* denotes the dual space of H^''^ . 
Then there exists a unique function v such that: 

(4.1) t;GL~(0,r;Fi'^), ^ e L\0,T; H^'^); 

dv \ , . 

— ,w-vj +ax{t;v,w -v)>0, 

\/w G H^'^ such that w > i/j{t) a.e. t£ (0, T); 

(4.3) v{■,t)>^p{t) a.e. te{0,T); 

(4.4) v{-,0) = v. 

Moreover, if v ^ L^(0, T; ff^'"^), then v is a solution to the obstacle problem: 
find V G L'^{0,T;H'^'^) such that v satisfies (4.3), (4.4) and 

(4.5) ^^+A{t)v>0; 

(4.6) (^ + A{t)v){v-i;) = 0, 



(4.2) 



.dt 

almost everywhere in M x (0, T) . 

Proof. For the most part, the theorem is a restatement of Bensoussan 
and Lions [1], Theorem 2.2, and Section 2.15, Chapter 3, where we have 
mapped t^T — t, and v 1— t- —v. 

We therefore only need to explain the last part of the result. If we suppose 
V G L^{0,T;H^'^) and cp G H^'^, we have 



ax{t;v,cP)= I e-2^Na(x,t)|^d<A+ / e'^^N^ 



bix,t)- 



dx 
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OV 



-OO 



+ / e-^^^'-'^<P-A{t)vdx, 



where the first term on the right-hand side vanishes since v G L^(0, t; i?^'^) 
and (j) £ H^'"^. Therefore, by (4.2), for any w G H^'''^ such that w>ip a.e. 
in R, 

dv \ 
— + A{t)v,w-v\ >0 a.e. t. 

/A 

Taking, for example, w = v + (j), ioi a, positive test function (/>, we conclude 
that (4.5) holds. Moreover, let tt; = in the inequality above, we have 

Then (4.6) follows from (4.3) and (4.5). □ 

4.2. Connection with Skorokhod's embedding problem. To connect our 
embedding problem SEP (u, i/, /x) with the variational inequahty, we need 
some assumptions on a, /j, and the starting distribution u. First, on aiM— t- 
M+, we still assume (3.2) and (3.4) hold. In addition, we assume that 

(4.7) 3K>0 such that ^ < a < on M. 

On /i and zv, we still assume that U^(x) < lJu{x) to ensure the existence of 
a solution to SEP (u, z/, //) . 

Under these assumptions, we can specify the coefficients in the evolution- 
ary variational inequality, (4.4) and (4.5)-(4.6), to be 

a{x, t) = ; b{x, t) = a{x)a'{x) — Xa'^{x) sgn(x); 

(4.8) 

V'(x,t) = U^(x); v = l]„{x), 
and then the corresponding operators are given by A{t) = and 



ax{t;v,w) 



-2X\x\ 



a'^(x) dvdw , \ 2r \ / ^^^^ 

+ [a[x)a [x) — Xa [x)sgi].{x))—w 



2 dx dx dx 
We write the evolutionary variational inequality as: 



dx. 



VI((T, V, ji): Find a function -y : M x [0, T] M satisfying (4.1)-(4.4), where 
ah the coefficients are given in (4.8). 

We also have a stronger formulation, that is: 

SVI(cr, i^, //): For given T > 0, we seek a function v, in a suitable space, 
such that (4.3)-(4.6) hold, where all the coefficients are given in (4.8). 
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Our main result is then to show that finding the solution to SEP(o", z^, ^) 
is equivalent to finding a (and hence the unique) solution to Yl{a,i', fi): 



Theorem 4.2. Suppose (3.2), (34) and (4.7) hold, and let T > 0. Also, 
let D and v be the solutions to SEP[a,v, fi) and VI{(t,v,^), respectively. 
Define u{x,t) := —W\x — XiatuI o.'^d by 

(4.9) D'^ :={(2;,t) gM X [0, T]; u(a;, t) > V(x, t)}. 

Then we have = D n M x [0, T] , and for all (x, t) G M x [0, T] , 

u{x, t) = v{x, t). 

Moreover, if u £ L^{0,T;H^''^), then u is also the solution to SVI{a,v,^). 



Proof. Let A > be fixed, and suppose D \s a solution to SEP (cr, u,^). 
We need to show u is a solution to W{a,v, First note that U^(2;) + 
is continuous on M, and converges to as x — >• ±00, and hence is bounded. 
Sox^V^{x) + \x\ gL~(0,T;/70'^), and then U^(x) G L~(0, T; i/O'^). Sim- 
ilarly, U^(x) G L°°(0,T;i?0'^). Since > U^(x) > u{x,t) > \J^{x) for all t G 
[0,r], we have u G L°°(0, T; ff°'^). By (3.6), we also have ||^| < 1 since 
u is the potential of some probability distribution. Therefore we have u G 
L°°{0,T; H^''^). By Lemma 3.3 and the fact that u is constant (in time) 
outside D, |^| < a'^p^{x,t) a.e. on M x [0,T] where p'^{x,t) is the transition 
density of the diffusion process X starting from i^. Then by standard Gaus- 
sian estimates (e.g., Stroock [36], Theorem 3.3.11), we know there exists 
some constant A>0, depending only on K, such that 



du 



dt 



L2{0,T;H0,^) 

-f ^ TT^'^^pl-^f^^-^^lT") -2\\x\\dxdtv{dy) 







T A ry+At 

e-2^l^ldxdti/(dy) 



lAtjj^_^j 

1^ A^ I ^^J_ '^i^-y? _2X\x\\dxdtu{dy) 

JrJo 1 * JM.\iy~At,y+At) [ 
T A rAt 



< 



[ [ [ e~^^^^Uxdtiy{dy) 



I r f f 2(x - y)^ , , , , , , 

+ / / — — / exp - ^ J^' \dxdtu{dy) 

J- /\ * JB\{y-At,y+At) I 
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1 



< 



1 At 
A f^2AXt 



(1 



-2XAt 



rT ^2At roo r 



2z^ 
' At 



dzdt 



1 At 



dt + 



^3/2^1/2 



^/2 



I At 



dt < oo, 



where we have apphed Holder's inequaUty in the first hne to get 



du 



dt 



p{t,y,x)u{dy) 



< 



p{t,y,xfu{dy). 



So If e L2(0,T;i/0'^), and we have shown (4.1) holds. 

By the same arguments used in the proof of Theorem 3.1, (4.3) and (4.4) 
hold. Now we consider (4.2). We begin by observing that, for any (j) G C^, 
if we write /X((dx) for the law of XtAr^i have 

d(t) du 



f d4>du ^ I 
dx dx 



dx 



(l-2P(XMr^ <x))dx 



(4.10) 



/ / ^'^{y<x}fJ-t{(iy)dx 
Jr Jr ox 



2E[(/>(XMr^)]. 



In addition, for any w G H ' , we can find a sequence {4>n} C such that 



(4.11) 



lim I 

n— ^oo 



{w-u{-,t) 



0. 



Moreover, e ^l^ln(x,t) is bounded, and if e ^\^\w is also bounded, then we 
can, in addition, find a sequence C such that e-2^l^l(/>„(x) > -K' 
for some constant K' independent of n. For any n, we therefore have 



(4.12) 



2 dx dx 



dx ■■ 



-2X\x\ 



{aa' - Acj^ sgn(x)) dx 



+ 



-2A|a;| 



c|)nO■'^^^Mx). 



On the other hand, since du/dt vanishes outside D, and, using the same 
arguments as (3.7) (which still hold on account of Lemma 3.3), is equal to 
-a{xfp^{x,t), we have, for almost every t G [0, T] 



(4.13) 



/ 

JR 



-2\\x\ 



dt 



-2X\x\ 



i„cr^/it(dx) 



-2X\x 



'Vn.O" /^((dx) 
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where Ds := {x eR: {x,s) € D}. By (4.12) and (4.13), 
'du 



dt 



+ ax{t;u,<pn) 

A 



--2A|a;| 



du a'^dud(j)n , , 2 / ^^^u 



'n 



JM. c/r Jig 

2A|x|i _2. 



/ 

Jr\d 



Utidx), 

for almost every t G [0,T]. Now suppose initially we have e~'^^^^w bounded, 
and choose a sequence as above. Then we can let n — >• 00 and apply 
Fatou's lemma and the fact that u = ij) on R \ and u; > -0 to get 

-(^,w-u\ + ax{t;u,w - u) = [ e'^^l'^l (w - ^)aVt(dx) > 0, 

V /a JM.\Dt 

for almost every t € [0,T]. So (4.2) holds when e~^^^^w is bounded. The 
general case follows from noting that maxji/;, — A^} converges to w in H^''^. 
We can conclude that u is a solution to VI(cr, i^, /i). In addition, the final 
statement of the theorem now follows from Theorem 4.1. 

Conversely, suppose that we have already found the solution to VI(cr, i^, //), 
denoted by v{x,t). By Theorem 4.1 and the preceding argument, we have 

-E'^lx - XtArol = V{x,t), 

when {x,t) € M x [0,T]. Finally, we need only note (from (3.7), and the 
line above) that whenever (x,t) G D, we have u{x,t) > ip(x,t), and hence 
D'^ = DnWx[0,T]. □ 

Remark 4.3. The constant A which appears in the variational inequal- 
ity can now be seen to be unimportant: if we consider two positive numbers 
A < A*, then by Theorem 4.1, there exist v and v* satisfying (4.1)-(4.4) with 
the parameters A and A*, respectively. According to Theorem 4.2, 

u{x,t) =v{x,t) =v*{x,t), 

so V = V* . Therefore, the description of Root's barrier by the strong varia- 
tional inequality is not affected by the choice of the parameter A > 0. We do, 
however, need A > 0, since this assumption is used, in, for example, (4.12), 
to ensure we can integrate by parts. 

Remark 4.4. As noted in Bensoussan and Lions [1], and which is well 
known, one can connect the solution to the variational inequality VI((T, z^, /x) 
to the solution of a particular optimal stopping problem. In our context, 
the function v which arises in the solution to VI(cj, z^, /x) is also the function 
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which arises from solving the problem 

(4.14) vix,t) = snpK'[\J^iXr)l{r<t}+^u{Xr)l{r=t}]- 

T<t 

This seems a rather interesting observation, and at one level extends a num- 
ber of connections known to exist between solutions to the Skorokhod em- 
bedding problem, and solutions to optimal stopping problems (e.g., Peskir [29], 
Obloj [28] and Cox, Hobson and Obloj [10]). 

What is rather interesting, and appears to differ from these other situ- 
ations, is that the above examples are all cases where the same stopping 
time is both a Skorokhod embedding, and a solution to the relevant optimal 
stopping problem. In the context here, we see that the optimal stopping 
problem is not solved by Root's stopping time. Rather, the problem given 
in (4.14) runs "backwards" in time: if we keep t fixed, then the solution to 
(4.14) is 

TD = inf{s >0:{Xs,t-s)^D}At. 

In addition, our connection between these two problems is only through 
the analytic statement of the problem: it would be interesting to have a 
probabilistic explanation for the correspondence. 

Remark 4.5. The above ideas also allow us to construct alternative 
embeddings which fail to be uniformly integrable. Consider using the varia- 
tional inequality to construct the domain D in the manner described above, 
but with the function •0 chosen to be V^{x) — a, for some a > 0. By (4.14), 
one can check that the solution to the variational inequality is a decreas- 
ing function with respect to t, and hence, B = is a barrier, which is 
nonempty, so that td < oo a.s., and the functions u{x,t) and v{x,t) defined 
in Theorem 4.2 agree (e.g., by taking bounded approximations to D). In 
particular, limt^oou{x,t) =\J^{x) — a. Since Xt/\rj^ is no longer uniformly 
integrable, we cannot simply infer that this holds in the limit, but we can 
consider for example 

U{X, t) - U{Z, t) = -E[\XtArD -X\- \XtATD " ^|] 

which is a bounded function. Taking the limit as t — )■ oo, we can deduce that 

-E[\Xr^ -x\- \Xr^ - z\] = U^(x) - U^(z). 

From this expression, we can divide through by {x — z) and take the limit 
as X 4, z to get 2F{Xrjy > z) — 1. The law of X^-^ now follows. 

Note also that there is no reason that the distribution above needed to 
have the same mean as z^, and this can lead to constructions where the 
means diff'er. In general, these constructions will not give rise to a uniformly 
integrable embedding, but if we take two general (integrable) distributions. 
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there is a natural choice, which is to find the smahest a G M such that 
U(/(a;) > U^(x) — a. In such a case, we conjecture that the resulting con- 
struction would be minimal in the sense that there is no other construction 
of a stopping time which embeds the same distribution, and is almost surely 
smaller. See Monroe [25] and Cox [8] for further details regarding minimality. 

4.3. Geometric Brownian motion. An important motivating example for 
our study is the financial application of Root's solution described in the 
Introduction. In both [16] and [5], the case a{x) = x plays a key role in 
both the pricing and the construction of a hedging portfolio. However, in the 
previous section, we only discussed the relation between Root's construction 
and variational inequalities under the assumptions (3.2), (3.4) and (4.7), 
where the last assumption is not satisfied by a in this special case. 

In this section, we study this special case: a{x) = x, so that Xt is a geomet- 
ric Brownian motion. In addition, we will assume that the process is strictly 
positive, so that u and /Li are supported on (0,cxd). We therefore consider 
the Skorokhod embedding problem SEP (cr, z^, /i) with starting distribution 
L', where v and n are integrable probability distributions satisfying 

supp(/i) C (0,oo), supp(i^) C (0,oo), 

(4.15) 



U^(2;) < Ujy(3;) and J x^ dv < 



CO. 



We recall from (2.1) that this implies, in particular, that the means of fx and 
V agree. 

The solution to the stochastic differential equation 

dXt = XtdWt, Xo = xo 

is the geometric Brownian motion XQexp{Wt — t/2}, and, for y > 0, the 
transition density of the process is 

. , 1 1 r (\nx-\ny + t/2f] 

By analogy with Theorem 3.1, if D is the solution to SEP((T, z^, /x), then 
we would expect 

at 9x2 u{x,t)=U^{x) on M X (0, oo) \ L»; 

where u is defined as before by u{x, t) = — E|x — Xiatd I • However, if we follow 
the arguments in Section 4.2, we find that we need to set a{x,t) = x'^/2 in 
VI((T, i^, /i), which would not satisfy the first condition of Theorem 4.1. To 
avoid this we will perform a simple transformation of the problem. We set 

v{x,t)=u{e'',t), (x,t) gM X [0,T]. 
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Define the operator A{t) := —\-§^ + \-§^] then we have, when {e^ ,t) G D, 

(4.17) ^^+A{t)v = ^. 

We state our main result of this section as follows: 



Theorem 4.6. Suppose a{x) = x on (0,oo) and /i and v satisfy (4-15). 
Moreover, assume D solves SEP{a,i', fx) , andu{x,t) := —M\x — XtArol- Then 
v{x,t) ■.= u{e^,t) is the unique solution to (4.1)-(4-4) where we set 

a(x,t) = i; 6(a;,t) = i-A-sgn(x); i) = U^(e^); 

(4.18) 

7; = U^(e^); A>i. 



Proof. Much of the proof will follow the proof of Theorem 4.2. As 
before, (4.3) and (4.4) are clear. In addition, we note that ip — e^ is continuous 
and converges to as x — )• oo and converges to U^(0) < oo as x — )• — oo, so 
x^ip-e"" £L°°{0,T;H°'^). Hence V S L°°(0, T; i7°'^) since we have A > ^ . 
Thus, V G L°°(0,r;ii'0'^). Moreover, we can easily see \dv/dx\ is bounded 
by e^. Therefore, v E L°°{0,T; H^'^) when A > ^. On the other hand, since 
\dv/dt\ is bounded by e^^ / pt(y, e^)i^(dy), we have, by Holder's inequality, 



dv 



dt 



< 



and hence, 
dv 



dt 



L2{0,T;H0-^) 
T 



< 



< 







T 



2TTt 



exp 



lny + t/2f 



+ 2x }v{dy), 



-2X\x\ 
2TTt 



■ exp 



(x-lny + t/2)^ 



27rt 



exp 



lny-t/2)^ 



+ 2x > dxdtiy{dy) 



t 



+ 2 In y > dx dt z^(dy) 



< 



y'^v{dy) 



1 



2J-Kt 



dt < oo. 



Therefore (4.1) is verified. 

Using (4.10), for (p £ we get 



(4.19) 



°° d du 

— [(j){ln{y))y] — {y,t)dy 
dy dx 

2E[cPilniXt,,rn))XtAro], 
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and so we define the measure vt by 

j <l){x)vt{dx)=^mn{Xtr,ru))XtArn]- 

Now take any w G ff^'-^, and take {(pn} C satisfying (4.11). By (4.17) 
and (4.19), similar arguments to those used in the proof of Theorem 4.2 give 

and 

f e-2A|x|^0^drE+ I e-^^\-\Mdx) = [ e-^^^^^M^x), 

JR (^t Jig JM.\Dt 

for almost all t G [0,r], where Dt := {x G M: (e"^,*) G D}. Thus, for almost 
every t G [0, T], 

f ( dv ^ 1 /I ^ / \A , f^^A 1 

= 1 U^" + 2 ^ 9^ + U " ''"^"^ J ^" 9^ J 

= /" e-2^l^l(^„i/i(dx). 

Finally, following the same arguments as in the proof of Theorem 4.2, we 
conclude (4.2) holds. Therefore f is a solution to (4.1)-(4.4) with coefficients 
determined by (4.18). The uniqueness is clear since it is easy to check the 
coefficients defined in (4.18) satisfy the conditions in Theorem 4.1. □ 

5. Optimality of Root's solution. For a given distribution /i, Rost [35] 
proves that Root's construction is optimal in the sense of "minimal residual 
expectation." It is easy to check that this is equivalent to the slightly more 
general problem 

minimize E[i^(r)] 
subject to: C{Xt) = /x; 

r is a UI stopping time. 

Here we assume /i is a given integrable and centered distribution, X is the 
diffusion process defined by (3.1), where the diffusion coefficient a satisfies 
(3.2)-(3.4), with initial distribution C{Xq) = v, and F is a given convex, 
increasing function with right derivative / and -F(O) = 0. 

Our aim in this section is twofold. First, since Rost's original proof re- 
lies heavily on notions from potential theory, to give a proof of this result 
using probabilistic techniques. Second, we shall be able to give a "pathwise 
inequality" which encodes the optimality in the sense that we can find a 
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submartingale Gt, and a function H[x) such that 

(5.1) F{t)>Gt + H{Xt) 

and such that, for To, equahty holds in (5.1) and Giato is a UI martingale. 
It then follows that td does indeed minimize EF(r) among all solutions 
to the Skorokhod embedding problem. The importance of (5.1) is that we 
can characterize the submartingale Gt, which will correspond in the finan- 
cial setting to a dynamic trading strategy for constructing a sub-replicating 
hedging strategy for call-type payoffs on variance options. 

We first define the key functions G{x,t) and H{x), where the submartin- 
gale in (5.1) is Gt = G{Xt,t), and give key results concerning these functions. 

We suppose that we have solved Root's problem for the given distribu- 
tions, and hence have our barrier B = D^. Define the function 

(5.2) M(x,t)=E(-'*)/(rB), 

where td is the corresponding Root stopping time. In the following, we shall 
assume 

(5.3) M{x, t) is locally bounded on M x M+. 

We suppose also (at least initially) that (3.2)-(3.4) and (4.7) still hold. Note 
that M(x,t) now has the following important properties. First, since / is 
right-continuous (it is the right derivative of F), M(x,t) = f{t) whenever 
(x, t) ^ D and t > 0. In addition, since / is increasing, for all x and t we 
have M{x,t) > f{t). 

Now define a function Z{x) by 

Jo Jo 

So in particular, we have Z"(x) = 2 ^I'^^'j^ , and Z{x) is a convex function. 
Define also 

(5.5) G{x,t)= [ M{x,s)ds- Z{x), 

Jo 

and 

rRix) 

(5.6) H{x)= / {f{s)-M{x,s))ds + Z{x), 

Jo 

where R{x) is the barrier function. Two key results concerning these func- 
tions are then: 

Proposition 5.1. We have, for all {x,t) G M x M+, 

(5.7) G{x,t) + H{x)<F{t). 



And also: 
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Lemma 5.2. Suppose that f is bounded, and for any T >0, 



(5.8) E 



/ Z'{Xsfa{X,fds 
Jo 



< oo, EZ{Xo) < oo. 



Then the process 

(5.9) G^XtATu^t Atd) is a martingale, 
and 

(5.10) G{Xt,t) is a submartingale. 

Using these results, we are able to prove the following theorem, which 
gives us Rost's result regarding the optimality of Root's construction. 

Theorem 5.3. Suppose D solves SEP(a, fi,iy), and equations (5.3) and 
(5.8) hold. Then 

(5.11) EF(rD)<EF(T) 
whenever t is a stopping time such that Xr ~ /-f. 

Proof. We begin by considering the case where Etd < oo,Et < oo and 
/ is bounded. Since Z{x) is convex, by the Meyer-Ito formula (e.g., Prot- 
ter [30], Theorem IV.71), 

Z{Xt) = Z{Xo) + / Z'{Xr) dXr + l [ Z"{Xr)<j\Xr) dr. 

Jo ^ Jo 

By (5.8) and the fact that / is bounded (and hence also M(Xc,, 0) is bounded), 
we get 

r-t/\T 

EZ(XtAr) =IE^(^o) +IE / M(X3,0)ds < /(oo)Er + EZ(Xo). 

Applying Fatou's lemma, we deduce that for any stopping time r with finite 
expectation, Z{Xr) is integrable. Moreover for such a stopping time, by 
convexity, Z{Xti\r) < E[Z(Xt-)|J^(], and so, by Lemma 5.2, G{Xt/\T,t ^t) is 
a submartingale which is bounded below by a UI martingale, and bounded 
above by /(oo)t. It follows that KG{XtAr,tAT) EG(X^,r) as t ^ oo. The 
same arguments hold when we replace r hy td. 

Since R{Xt-^) < td and if t € [i?(x),oo), then td = t,P(^''*^-a.s., so that 
M{Xt-jj , s) = f{s) for s>T£), we have 

G{Xr^ ,Td)+ ifis) - M{Xr^ , S)) ds + Z{XrJ 
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/ 

Jo 



L 



M{Xra,s)ds + 



{f{s)-M{X^^,s))ds 



(5.12) 




/ M{Xr^,s)ds+ {f{s)-M{Xr^,s))ds 



Jo Jo 



/ f{s)ds = F{TD). 



Jo 



On the other hand, since X^-^ ^ Xr, and observing that G{Xrjy,TD) and 
F{t]:)) are integrable, so too is H{Xt-jj), and 



In addition, by Lemma 5.2 and the hmiting behavior deduced above, we 
have 



We now consider the case where at least one of r or td has infinite expec- 
tation. Note that if F{-) ^ 0, then there is some a, /3 G M with /? > such that 
F{t) >a + f3t, and hence we cannot have Er = oo or Er/j = oo without the 
corresponding term in (5.11) also being infinite. The only case which need 
concern us is the case where Er < oo, but Ktd = oo. Note, however, that td 
remains UI, so E[XtAT-D |-^t] = Xt. In addition, from the arguments applied 
above, we know ZiX^) is integrable, and since X^- ~ X^-^^^ so too is Z(Xt-^). 
Then H{Xt-) and H{Xrj-i) are both bounded above by an integrable random 
variable, so their expectations are well defined (although possibly not finite), 
and equal. Then, as above, — E[Z(Xt-^)| 7"^] < —Z{XtATu) ^ G{Xt/\Tjj,tATD)- 
We can deduce that EG(X^^ ,td) < lim„^oo EG(XtAro ,t Atd) = G{Xo,0) < 
KG{Xt-,t). The remaining steps follow as previously, and it must follow that 
in fact EF(r/)) < EF(r), which contradicts the assumption that Er < oo and 
Ktd = oo. 

To observe that the result still holds when / is unbounded, observe that 
we can apply the above argument to f{t) A A^, and F^it) = Jq f{s) A N ds 
to get EF/v(tz)) < EF/v(t), and the conclusion follows on letting N ^ oo. 
□ 

We now turn to the proofs of our key results: 



EH{XrJ=EH{Xr). 



EG{Xro,TD) = EG{Xo,0) < lim EG{XtAr,t At) = EG{Xr,T). 



Putting these together, we get 

EF{td) = E[G{Xr^,TD) + H{Xr^)] 
<E[G{Xr,T) + H{Xr)] 

<EF{t). 
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Proof of Proposition 5.1. If i < R{x), then the left-hand side of 
(5.7) is 

ff{s)ds+ r^^\f{s)-M{x,s))ds = F{t)- r^^\M{x,s)-f{s))ds, 
Jo Jt Jt 

and we know M{x, s) > f{s) > 0, so that the inequahty holds. 

Now consider the case where R{x) < t. Then the left-hand side of (5.7) 
becomes 

i-R(x) pt pR(x) 



/ M{x,s)ds+ / f{s)ds= / f{s)ds+ / f{s)ds = F{t). 
Jr(x) Jo Jr(x) Jo 



□ 



Proof of Lemma 5.2. We begin by noting that Z{x) is convex, and 
therefore the Meyer-Ito formula (e.g., Protter [30], Theorem IV. 71) gives 

Z{Xt) - Z{Xs) = Z'{Xr) dXr + ^J^ Z"{Xr)a\Xr) dr. 
It follows from (5.8) that the first integral is a martingale. So we get 

E[Z{Xt) - ZiX^m = [ E[M{Xr,0)\Ts]dr, s < t. 

J s 

In addition, since M{x,t) > f{t) and /(t) is increasing, for r, u > by the 
strong Markov property, writing X for an independent stochastic process 
with the same law as X and td for the corresponding hitting time of the 
barrier, we have 

]E("'^)[/(r,,)|J-,+J = lr,>r+uE^-^'\f{rD)\Fr+u] 

+ l,^<,,+„E(^''^)[/(rz))|j;+„] 

< l,^>,+„E(^S'^+«)[/(fB)] + lru<r+uf{r + U) 

<M{X:,r + u). 

When r = 0, we have E^^'O) [/(rD)|-F«] <M{X^,u). For s,nG [0,t], 
K[M{Xt,u)\Ts]=K^'M{Xt_s,u) 

(5.13) >e(^-"-(*-"))[/(?d)] 

>A'I{Xs,u-{t-s)), 

when u>t — s. On the other hand, if n < t — s, 

E[M{Xt,u)\Ts] = E[e(^'— 0) [M{Xu,u)]\Ts] 

(5.14) >E[E(^'-"'0)[/(fB)]|-F,] 

>E[M(Xt_„,0)|J-,]. 
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Then we can write 

ft 



E[G{Xt,t)\Ts]= / nMiXt,u)\J^s]du-nZ{Xt)\J^s] 
Jo 

= G{Xs,s)+ [ E[M{Xt,u)\J^s]du- [ M{Xs,u)du 
Jo Jo 

-E[Z{Xt)-Z{Xs)\Ts] 

>G{Xs,s)+ I nM{Xt-uM^s]du- r M{Xs,u)du 
Jo Jo 

- I E[M{XuM^s\du+ j M{Xs,s-t + u)du 

J s Jt—s 

>G{Xs,s)+ fE[M{Xu,0)\Ts]du- fE[M{Xu,0)\Ts]du 

J s J s 

+ / M{Xs,u)du- / M{Xs,u)du 
Jo Jo 

>G{Xs,s). 

Where we have used (5.13) and (5.14) in the third hne. 

On the other hand, on {td > s}, from the definition of M{x,t) and the 
Markov property, we get 

(5.15) E[M{Xt^ro,tA TD - u)\Ts] = M{Xs,s - u) 
when u< s, and 

(5.16) E[M(XMr« ,tATD- n)| = M(X„, 0) 

when u£ [s,t A td]. Then a similar calculation to above gives, for s <to-, 
E[G{Xt/,ru.t^rD)\Fs] 



■-E 



thru 



M{XtATD ,t Atd - u)du 







E[ZiXt;,rn)\J'si 



M{Xs,s-u) du + E 

-tATD 



t/\TD 



M{XtATD,'t^TD-u) du 



- Z{X,) - E 

ft 



M{Xu,0) du 



■-E 



I E[M {Xt^ro ,tATD-u)- M{Xu,0)\J'u]l{u<ro} 
.J s 

+G{Xs,s) 
= GiX,,s), 
where we have used (5.15) and (5.16). □ 
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Remark 5.4. Note that the fact that our choice of D given in the 
solution is the domain D which arises in solving Root's embedding problem 
is only used in Theorem 5.3 to enforce the lower bound. In fact, we could 
choose any barrier B, and D = as our domain, and this would result in a 
lower bound, with corresponding functions G and H. The choice of Root's 
barrier gives the optimal lower bound, in that we can attain equality for 
some stopping time. In this context, it is worth recalling the lower bounds 
given by Carr and Lee [5], Proposition 3.1 — here a lower bound is given 
which essentially corresponds to choosing the domain with R{x) = Q, for a 
constant Q. The arguments given above show that similar constructions are 
available for any choice of R, and the optimal choice corresponds to Root's 
construction. 

Remark 5.5. Although the preceding section is written for a diffusion 
on M, it is not hard to check that the case where cr{x) = x can also be included 
without many changes. In this setting, we need to restrict the space variable 
to the space (0,oo) (so we assume that r^i < oo a.s.), and consider a starting 
distribution which is also supported on (0,oo), and with a corresponding 
change to (5.3). 

We end this section with a brief example which illustrates some of the 
relevant quantities. 



Example 5.6. Suppose we take Root's barrier D := {{x,t):t < R{x)} 
with the boundary function R{x) = — A(j; + q;)(x — /3)1(_q,j3), where A,a,/3 > 
0; see Figure la. Given a standard Brownian motion W and Root's stopping 
time TD = inf{t > 0:t > RiWt)}, define n := C{Wro)- Let F{t) = and 
we will see E[F(r£))] < E[F(r)] for any UI stopping time r such that Wr ~ /i. 

For (x, t) G R X M+, define M{x,t)= E^^-*) [td] ■ Then iit> R{x), M{x,t) = 
t. If < t < R{x), since td = A(WVr, + ce)(l^TD ~ f^)^ using Ito's formula, we 
can compute M{x,t) to be 

M{x,t) = -^—[t - {x + a){x - f3)] for 0<t<R{x). 
1 + A 

Defining G,H,Z as in (5.4)-(5.6), we get the explicit expressions 



Z{x) 



G{x,t) 



A 



6(1 + A) 



2a/33 + (2/33 + 6a/?2)2;, 
x^-2(a-/3)x3 + 6a/3a;2. 



-a 



A 



1 + A 

R^{x 



t{x + a){x — /3) 



Z(x), 



1 2 

2(1 + A) ^2* 



Z(x), 



x>/3, 

X G (— a,/3), 

X < —a, 

[fO<t<R{x), 
if t>R{x), 
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Fig. 2. We give graphical representations of some of the relevant quantities derived in 
Example 5.6, for a = 2, /? = 3 and A = 1/2. In (a) we see G{x,t) + H{x), which is a lower 
bound for F{t), and in (b) we see the difference G{x,t) + H{x) — F{t), which is indeed 
negative. 



H{x) 



R\x) 



■Z{x). 



2(1 + A) 

It is easy to check directly that G{Wt,t) is a submartingale, and that it is a 
martingale up to the stopping time td- We also can check that (5.7) holds 
here: 

[R{x)-tf 



G{x,t)+H{x)-F{t) 



2(1 + A) 



if 0<t <i?(x), 

0, ift>R{x). 
Therefore, for any UI stopping time r such that C{Wr) = ^ = ^(W^r^)) 

E[F(r)] > E[G{Wr,T) + H{Wr)] > E[G{Wr^,TD)] +E[H{Wr^)] 

(5.17) 



:E[F{R{Wr^))]+K 



M{Wro,s)ds 



nFiro)], 



which shows the optimality of Root's stopping time. Figure 2 illustrates 
some of the relevant functions derived here. 



6. Financial applications. We now turn to our motivating financial prob- 
lem: consider an asset price St defined on a complete probability space 

iTt)t>o,w), with 

(6.1) ^ = rtdt + atdWt 

^t 

under some probability measure Q ~ P, where P is the objective probability 
measure, and Wt a Q-Brownian motion. In addition, we suppose rt is the 
risk-free rate which we require to be known, but which need not be constant. 
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In particular, let rj,(Tj be locally bounded, predictable processes so that the 
integral in (6.1) is well defined, and so St is an ltd process. We suppose 
that the process at is not known (or more specifically, we aim to produce 
conclusions which hold for all at in the class described) . Specifically, we shall 
suppose: 

Assumption 6.1. The asset price process, under some probability mea- 
sure Q ~ P, is the solution to the SDE (6.1), where rt and at are locally 
bounded, predictable processes. 

In addition, we need to make the following assumptions regarding the set 
of call options, which are initially traded: 

Assumption 6.2. We suppose that call options with maturity T, and 
at all strikes {K : K > 0} are traded at time 0, and the prices, C{K), are 
assumed to be known. In addition, we suppose call-put parity holds, so that 
the price of a put option with strike K is P{K) = "^sdsj^ — 5o -|- C{K). 
We make the additional assumptions that C(K) is a continuous, decreasing 

and convex function, with C(0) = Sq, C'^{0) = -e~fo'-^'^' and C{K) ^ as 
K ^oo. 

Many of these notions can be motivated by arbitrage concerns; see, for 
example, Cox and Obloj [12]. That there are plausible situations in which 
these assumptions do not hold can be seen by considering models with bub- 
bles (e.g., [9]), in which call-put parity fails, and C{K) 7^ as — >■ oo. Let 

us define Bt = e^o and make the assumptions above. Following the per- 
spective that the prices correspond to expectations under Q, the implied 
law of B^^St (which we will denote /i) can be recovered by the Breeden- 
Litzenberger formula [2], 

(6.2) ^i{{K, 00)) = q*{B~^ St e {K, 00)) = -2BtC'_^{BtK). 

Here we have used Q* to emphasize the fact that this is only an implied 
probability, and not necessarily the distribution under the actual measure 
Q. From (6.2) we deduce that U^(a;) = — 2C{Btx) — x, giving an affine 
mapping between the function Uo(x) and the call prices. We do not impose 
the condition that the law of B^ St under Q is fi, we merely note that this 
is the law implied by the traded options. We also do not assume anything 
about the price paths of the call options: our only assumptions are their 
initial prices, and that they return the usual payoff at maturity. It can 
now also be seen that the assumption that C^(0) = —B^^ is equivalent to 
assuming that there is no atom at — that is, /i is supported on (0,oo). 
Finally, it follows from the assumptions that /i is an integrable measure 
with mean Sq. 
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Our goal is to now to use the knowledge of the call prices to find a lower 
bound on the price of an option which has payoff 

FQ\Ut^=F{{lnS)T). 
Consider the discounted stock price, 

Xt = e-^o^^'^'St = B^^St. 
Under Assumption 6.1, Xt satisfies the SDE 

dXt = XtatdWt. 

Defining a time change Tt = /q <Jg ds, and writing At for the right-continuous 
inverse, so that = t, we note that Wt = ds dWs is a Brownian motion 
with respect to the filtration F't = ^At > c^nd if we set Xi = Xa_^ , we have 

dXt = XtdWt. 

In particular, Xt is now of a form where we may apply our earlier results, 
using the target distribution arising from (6.2), and noting also that Xq = Sq 

and Xrj, = Xt = Bj^^St- 

We now define functions as in Section 5, so that f{t) = F^{t) and (5.2)- 
(5.6) hold. Our aim is to use (5.7), which now reads 

(6.3) G{XAt , t) + H{XAt ) = GiXt ,t) + H{Xt)< F{t) = F a] d^ , 

to construct a sub-replicating portfolio. We shall first^show that we can 
construct a trading strategy that sub-replicates the G{Xt,t) portion of the 
portfolio. Then we argue that we are able, using a portfolio of calls, puts, 
cash and the underlying, to replicate the payoff H{Xt)- 

Since G{Xt,t) is a submartingale, we do not expect to be able to replicate 
this in a completely self-financing manner. However, by the Doob-Meyer 
decomposition theorem, andj^he martingale representation theorem, we can 
certainly find some process (j)t such that 

G{Xt,t)>G{Xo,0)+ [ ^sdXs 
Jo 

and such that there is equality at t = tq- Moreover, since G{XT-j^^f, td A t) 
is a martingale, and G is C^'^ in D, we have 

G(X^^At,TD At) = G(Xo,0) + —iXrj,As,TDAs)dXs. 

More generally, we would not expect ^ to exist everywhere in D^; however, 
if, for example, left and right derivatives exist, then we could choose (pt G 
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^(^"1"'*)] holding of the risky asset (or alternatively, but 

less explicitly, take = (9/9x[E^'*G(Xf+5, to + <5)], for t G [tQ,to + 6]). 
It follows that we can identify a process (pt with 

G{X^, ,Tt)> G{Xo, 0) + r dXs = G{Xo, 0) + /* dX^, 

Jo Jo 

where we have used, for example, Revuz and Yor [31], Proposition V.1.4. 

Finally, writing (pg = 4>ts > we have 

G{Xt,Tt)>G{Xo,0)+ f (l)sdXs = G{Xo,0)+ f c^sdiB^^Ss). 

Jo Jo 

If we consider the self-financing portfolio which consists of holding cpsB^^^ 

units of the risky asset, and an initial investment of G{Xq,0)B^^ — (poSoBj^^ 

in the risk- free asset, this has value Vt at time t, where 

d{Br'Vt) = B^'ct>td{Br'St), 

and therefore 

Vt = Bt (voB^' + £ 5^ V. diB^^Ss)^ = G{Xo,0) + £ (t>s dX^. 

We now turn to the H{Xt) component in (6.3). If H{x) can be written 
as the difference of two convex functions (so, in particular, H"[dK) is a 
well-defined signed measure), we can write 

H{x) = H{So) + H'4So){x -So)+ [ {x - K)+H"{dK) 

J{So,oo) 

+ [ iK-x)+H"{dK). 

J{0,So] 

Taking x = Xt = B^^St, we get 

H{Xt) = H{So) + H'_,{So){B^'St - So) + B^' [ {St - BTK)+H"{dK) 

J{So,oc) 

+ B^^ [ {BtK - ST)+H"{dK). 

J{0,So] 

This implies that the payoff H{Xt) can be replicated at time T by "holding" 
a portfolio of 

B~^{H{So) - H'^{So)So) in cash; 
B^^ H',{So) units of the asset; 

(6.4) 

B^'^H"{dK) units of the call with strike BtK for K G (5o,oo); 
B;^^H"{dK) units of the put with strike BtK for K £ (0,5o]; 
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where the final two terms should be interpreted appropriately. In practice, 
the function H{-) can typically be approximated by a piecewise linear func- 
tion, where the "kinks" in the function correspond to traded strikes of calls 
or puts, in which case the number of units of each option to hold is deter- 
mined by the change in the gradient at the relevant strike. The initial cost 
of setting up such a portfolio is well defined, provided 



where |-ff"|(d-ftr) is the total variation of the signed measure H"(dK). We 
therefore shall make the following assumption: 

Assumption 6.3. The payoff H{Xt) can be replicated using a suitable 
portfolio of call and put options, cash and the underlying, with a finite price 
at time 0. 

We can therefore combine these to get the following theorem: 

Theorem 6.4. Suppose that Assumptions 6.1, 6.2 and 6.3 hold, and 
suppose F{-) is a convex, increasing function with -F(O) = and right deriva- 
tive f[t) = F^{t) which is bounded. Then there exists an arbitrage if the price 
of an option with payoff F {{In S)t) is less than 



where the functions G and H are as defined in (5.5) and (5.6), and are 
determined by the solution tq to SEP{a,6so, ^j) for cr{x) =x, and where // 
is determined by (6.2). 

Moreover, this bound is optimal in the sense that there exists a model 
which is free of arbitrage, under which the bound can be attained. 

Proof. It follows from Theorem 4.6 that, given /x, we can find a do- 
main D and corresponding stopping time td which solves SEP((7, 5so,n). Ap- 
plying Proposition 5.1 (and bearing in mind Remark 5.5), we conclude that 
the strategy described above will indeed sub-replicate, and we can therefore 
produce an arbitrage by purchasing the option, and selling short the portfo- 
lio of calls, puts and the underlying given in (6.4), and in addition, holding 
the dynamic portfolio with —(f>tB^^ units of the underlying at time t. It is 
not hard to check, given that / is bounded (and choosing the lower limits 

in (5.4) to be 5o rather than 0) that {Z'{Xs)a{Xs)f < {Xs/Xq - l)^ and 
hence that (5.8) holds. Condition (5.3) also clearly holds. As a consequence, 
we do indeed have a subhedge. 



(6.5) 




(6.6) 
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To see that this is the best possible bound, we need to show that there is 
a model which satisfies Assumption 6.1, has law /x under Q at time T, and 
such that the subhedge is actually a hedge. But consider the stopping time 
td for the process Xt- Define the process 

Xt = Xt/{T-t)^ro fortG[0,r] 

which corresponds to the choice of al = ^3^1{s/(t-s)<td}- Since T£, <oo 

a.s., then Xt = X^-jj, tt = td and St = XtBf is a price process satisfying 
Assumption 6.1 with 

FQ\^dt^ =F{td). 

Finally, it follows from (5.12) that at time T, the value of the hedging 
portfolio exactly equals the payoff of the option. □ 

Remark 6.5. The above results are given in the context of an increas- 
ing, convex function, but there is also a similar result concerning increasing, 
concave functions which can be derived. Consider a bounded, increasing 
function / as before, and define the function 

L{t) = A/(oo) - f{s))ds = /(oo)t - F{t). 
Jo 

Using Theorem 6.4 and (1.2), it is easy to see that the price of a contract 
with payoff L{{\nS)T) must be bounded above by 

2/(oo)Q - 2/(oo)5^i log(5o) - ^^^^^(^o, 0) - 5^^i7(5o) 

-B^^ [ C{BTK)H"{dK)-B^^ f P{BTK)H"{dK), 

where Q is the price of a log-contract [i.e., an option with payoff ln(5T)]. 
As before, this upper bound is the best possible, under a similar set of 
assumptions. 

Remark 6.6. An analogous result can be shown for forward start op- 
tions. Suppose that the option has payoff 




= F{{S)t-{S)s) 



for fixed times < S <T. Then we can use the previous results for general 
starting distributions to deduce a similar result to Theorem 6.4 for forward 
start options, provided we assume that there are calls traded at both S and 
T. We use essentially the same idea as above: we aim to hold a portfolio 



ROOT'S BARRIER: CONSTRUCTION, OPTIMALITY AND APPLICATIONS 33 



which (sub-)rephcates G{Xt,Tt) for t G [S, T], and hold the payoff H{Xt) 
as a portfoho of calls. However, we now have rj = J^c^ds, and so Xt = 

XAt 1 gives Xq = Xs (recall that At was assumed right-continuous) . The 
procedure is much as above, except that we need to use the solution to 
Theorem 5.3 with a general target distribution, and the amount G(Xo,0) 
will be a J^5-random variable. The initial distribution v can be derived 
using the Breeden-Litzenberger formula (6.2) at time S. To ensure that we 
hold the amount G(Xo,0) at time 5, we observe that G(Xo,0) = G{Xs,Q)- 
Hence if, for example, G(x,0) can be written as the difference of two convex 
functions, we can replicate this amount by holding a portfolio of calls and 
puts with maturity S in a similar manner to (6.4). The remaining details 
follow as in the hedge described in Theorem 6.4 

Remark 6.7. We can also consider modifications to the realized vari- 
ance. Consider a slightly different time-change: suppose we set 

Tt= f al\{Xs)ds, 
Jo 

for some "nice" function X{x), which in particular we suppose is bounded 
above and below by positive constants. Then following the computations 
above, we see that 

Xt = XA, = r\sKXs)-^'\cjA{X,f^dWs)= f XaXXaX^'^AWs, 
Jo Jo 

and therefore dXt = cr{Xt) dWt, where a{x) = xA(x)~^/^. We then conjecture 
that it is possible to extend Theorem 4.6 to cover this new class of functions 
cr{x) (the conditions that should be imposed on A such that this result may 
be extended remains an interesting question for future research). It would 
then be possible to modify the above arguments to provide robust hedges 
on convex payoffs of the form 

F(^l\',XiXs)dsy 

An interesting special case of this would then be to give robust bounds on 
the price of an option on corridor variance 

(6.7) F(^l\ll{Sse[a,b]}<^^y 

by considering A(x) = l{a;g[a,6]}5 however this would only work in the case 
where there are no discount rates (i.e., Bt = l). In general, we can only give 
a tight lower bound for options on 
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although this does provide a lower bound for (6.7) by considering the case 
where a = a and b = Bxb. 

7. Conclusions. We conclude by summarizing the results, and describ- 
ing some interesting questions for future work. In this paper, we have given 
a variational inequality representation of Root's solution to the Skorokhod 
embedding problem, and provided a novel proof of optimality, which allows 
us to construct a model-independent subhedge for options on variance. We 
believe that our results provide interesting insights into all three aspects of 
the work: the construction of solutions to the Skorokhod embedding prob- 
lem, proving optimality results for the same and finally the connections with 
model-independent hedging. 

We also believe that there are interesting lines of research that now arise. 
The construction opens up a number of questions regarding Root's solution 
to the Skorokhod embedding problem: for example, what can be said about 
the shape of the boundary? Under what conditions on fi will the boundary 
be smooth? When does R{x) — )• as x — )• ±oo? When is R{x) bounded? 
Properties of free boundaries are well studied in the analytic literature, and 
may be useful in answering these questions. The connection to minimality 
and noncentered target distributions raised in Remark 4.5, and the question 
asked at the end of this remark would also be interesting lines for research. 

The connection with optimal stopping noted in Remark 4.4 is interesting, 
and obtaining a deeper understanding between optimal stopping problems 
and optimal Skorokhod embeddings seems to be an interesting area of re- 
search. 

Another natural question concerns the upper bound/super-hedging strat- 
egy. It has been remarked by Obloj [27] and Carr and Lee [5] that a related 
construction of Rost should provide a suitable upper bound, but similar 
questions to those answered here remain (although we hope to be able to 
provide some answers in subsequent work). We note, however, that numerical 
evidence (see Carr and Lee [5]) seems to suggest that the Root bounds may 
be more appropriate in the financial applications. It would also be of inter- 
est to see to what extent these model-independent bounds may be useful in 
practice. In Cox and Obloj [12], an analysis of the use of model-independent 
bounds as a hedging strategy for barrier options was performed. A similar 
analysis of the strategies derived in this work would also be of interest. 

Other questions that arise from the practical standpoint include how to 
incorporate additional market information (e.g., calls at an intermediate 
time [4]), and how to adjust for the fact that there will generally only be 
a finite set of quoted calls; see [13] for a related question. Remark 6.7 also 
suggests open questions regarding more general choices of a{x). 

Acknowledgment. We are grateful to Sam Howison for a helpful discus- 
sion which has much improved the material in Sections 3 and 4. 
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